Structure and Magnetic Fields in the Precessing Jet System SS 433 
III. Evolution of the Intrinsic Brightness of the Jets from a Deep 

Multi-Epoch VLA Campaign 



Michael R. BelV David H. Roberts, and John F. C. Wardle 
Department of Physics, MS-057, Brandeis University 
415 South Street, Waltham, MA 02453 
|mrbell@mpa-gar ching . mpg . de| , IrobertsQbrandeis . edul , iwardleQbrandeis . eduj 



ABSTRACT 

We present a sequence of five deep observations of SS 433 made over the summer of 2007 using 
the VLA in the A configuration at 5 and 8 GHz. In this paper we study the brightness profiles of 
the jets and their time evolution. We also examine the spectral index distribution in the source. 
We find (as previously reported from the analysis of a single earlier image) that the profiles of 
the east and west jets are remarkably similar if projection and Doppler beaming are taken into 
account. The sequence of five images allows us to disentangle the evolution of brightness of 
individual pieces of jet from the variations of jet power originating at the core. We find that 
the brightness of each piece of the jet fades as an exponential function of age (or distance from 
the core), e~^l^ , where r is the age at emission and r' = 55.9 ±1.7 days. This evolutionary 
model describes both the east and west jets equally well. There is also significant variation (by 
a factor of at least five) in jet power with birth epoch, with the east and west jets varying in 
synchrony. The lack of deceleration between the scale of the optical Balmer line emission (10^^ 
cm) and that of the radio emission (10^^ cm) requires that the jet material is much denser than 
its surroundings. We find that the density ratio must exceed 300:1. 

Subject headings: binaries: close — radio continuum: stars — stars: individual (SS433) 



1. Introduction 

SS 433 was the first known microquasar, a type 
of x-ray binary (XRB) system consisting of a col- 
lapsed star accreting material from a less evolved 
donor star. The feature that distinguishes a micro- 
quasar from other XRBs is that they are variable 
radio sources that can eje ct pairs of oppositely 
directed relativistic jets (Mir abel fc Rodrigued 
fl999l In many ways they are analogous to their 
AGN counterparts only smaller in size resulting in 
much shorter variability time scales. This provides 
a distinct observational advantage for studying ac- 
creting black holes and relativistic jet systems. 



SS433 has been intensely studied since the 
discovery of oscillating simultaneously red and 
blue s hifted hydrogen lines in the optical spec- 
trum (Margon et al.'"l979a'). A kinematic model 
was proposed that explained these features as 
the resul t of a pair of mildly- relativistic precess- 
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ing jets (Margon et ajj Il979bl ) . This model was 
soon confirmed and several parameters disam- 
bi^ua ted by comparison with early VLA radio im- 
ages ( Hjellmi ng fc J ohnstonI [l981[) . The jets pre- 
cess with a period of 163 days about a cone with 
a semi-opening angle of ~ 20°. The central axis is 
inclined from the line of sight by ^80°. For in - 
dep th reviews of pr evious work, see lMargonl ( 1984J ) 

and lFabrikal (|2004]). 

Previously in Roberts et al. ( 20081 ) we analyzed 
the structure of SS 433 in both total and polarized 
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intensity a t 0.1^^ resolution from a 15 GHz VLA 
image. In [Roberts et al.l ( 2Q1Q| ) (henceforth Pa- 
per II) we described in detail the procedure for 
measuring the intrinsic brightness of the jets, a 
measure of the source emission as observed in a 
comoving reference frame, and analyzed the pro- 
file of a single, high dynamic range observation. 
We found that the east and west jet profiles were 
consistent with an assumption of intrinsic symme- 
try and that the profile was more complex than 
could be fit by a simple exponential or power law. 
We suggested that the profile shape may be in- 
dicative of variation in the power injected into the 
jet. 

In this paper we present a series of deep two- 
frequency VLA observations of the jets of SS 433 
and use them to study the jets as they evolve over 
an 80 day period. Proper characterization of the 
jet dynamics is important for testing and devel- 
oping jet models which predict, amongst other 
things, the brightness decay rate based on as- 
sumed physical parameters. 

In Section [2] we describe the observations and 
data reduction and present the images. In Sec- 
tion [3] we present maps of the spectral index across 
the source. In Section |4] we calculate the intrinsic 
brightness profile at each epoch and investigate 
how it evolves in time. We discuss the implica- 
tions of our results in Section [5] and present our 
conclusions in Section [6l 

2. Observations 

In order to study the evolution of the arcsecond 
scale radio emission from SS433, we conducted a 
series of five VLA A-array observations made over 
the summer of 20070 Each observation includes 
data at two frequencies: 5 GHz (A6 cm, C-band) 
and 8.5 GHz (A3. 6 cm, X-band). Table [T] provides 
a summary of the observing dates and the preces- 
sional phases for the observations. 

Frequent observations of the nearby, unresolved 
calibrator source J 1950+081 were included for 
phase calibration purposes. Two fiux calibrator 
sources, 3C286 and 3C48, were observed at the be- 
ginning and end of the observations, respectively. 
Additionally, a 5 minute scan of the bright source 
3C84 (-14 Jy at 5 GHz as of 2007) was included 



^VLA project code AR0637 



to allow for the calibration of baseline-dependent 
complex gains. 

These observations were deemed "shared-risk" 
because the VLA was in the midst of being up- 
graded to the EVLA. This meant that the ar- 
ray was in a transitional state, with ten anten- 
nas having received new hardware while the oth- 
ers were still equipped with legacy hardware. The 
upgraded antennas had been equipped with new 
electronics, including the IF/LO system, receivers, 
and bandpass filters. Also, on June 27th, the 
"MODCOMP" computer system that was respon- 
sible for controlling the array was upgraded to the 
EVLA monitor and control system. The WIDAR 
correlator was not yet in use at this time. The 
number of VLA and EVLA antennas used in each 
observation, after data fiagging, is shown in Table 

m 

All data were calibrated using the NRAO As- 
tronomical Image Processing Software (AIPS) 
package. We followed the sta ndard procedu re out- 
lined in the AIPS Cookbook (InRAOI|20oI ). mod- 
ified to include baseline calibration. The baseline 
calibration procedure is not typically necessary 
when processing VLA data, but due to the mis- 
matched bandpass filters between VLA and EVLA 
antennas, this step was essential to reduce closure 
errors. Imaging and self-calibration were carried 
out using DIFMAP ([Shepherd et al.lll994l ). 

Due to last minute re-scheduling and problems 
with the array, the July 1st data at both frequen- 
cies and the August 24th data at 5 GHz did not 
include observations of a fiux calibrator. Instead, 
the fiux of the phase calibrator was determined 
from the other observations and used to set the 
absolute fiux scale. The fiux of J 1950+081 was 
roughly constant throughout the summer, measur- 
ing 1.09 ± 0.02 Jy at C-band and 0.78 ± 0.02 Jy at 
X-band. We estimate the uncertainty in our fiux 
scales as ±3%. 

2.1. The Images 

Figures [T] and [2] show the uniformly weighted 
total intensity images at each epoch for 5 and 8.5 
GHz, respectively. The plots are of total inten- 
sity with a step size of a/2 between contours. The 
CLEAN beam is shown as a cross in the lower right 
hand corner. More information about each image 
can be found in Tables [2] and [H In these tables the 
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columns marked "BMAJ" and "BMIN" provide the 
FWHM angular dimensions of the major and mi- 
nor axes of the restoring beam; the beam position 
angle is listed in the column marked "BPA". The 
Imax and Imin valucs are the map peak and bot- 
tom contour intensities, respectively. The "RMS" 
value is the average of the root mean squared in- 
tensity measured in four empty regions of the im- 
age. For some images the RMS has also been cal- 
culated by fitting a histogram of the pixel inten- 
sity values to a Gaussian. Both methods yield the 
same results. The thermal noise limit of the VLA 
for these observations as calculated using the VLA 
exposure calculatoiH is approximately 20/iJy for 
the C-band images and 14/iJy for the X-band im- 
ages. The RMS noise level achieved in each image 
is typically 2-3 times these values. 

The simple kinematic model is overlaid on 
each image with blue and red line segments in- 
dicating material that is traveling toward and 
away from the obs erver, respectively . We adopt 

(|2QQl[ ) and 
of 5.5 kpc 

(jBlundell fc Bowlerl 12004 ). A summary of the 
model parameters may be found in Table 1 
of Paper II. We include neither the nutation 



the ephemeris of lEikenberrv et al 
assume a distance to the source 



Katz et al.l Il982r) no r orbital velocity variation 



Blundeh etliD l2007h ince they imply structure 
on a smaller scale than the ^ 350 mas resolution 
of our images. In Paper II these perturbations 
to the basic model were included as a method of 
sampling the area around the mean model curve 
in order to estimate the degree to which our mea- 
surements varied by changing the model location. 
Here we use an alternate method for determining 
this uncertainty as described below. 

The core flux in each observation is shown as 
a function of truncated Julian date (TJD = JD - 
2454000) in Figure [3l The error bars in the figure 
represent the 3% calibration uncertainty. This is 
much larger than the uncertainty due to noise in 
the image and is a conservative estimate based on 
the stability of the measured flux of the phase cal- 
ibrator. We find that the 5 and 8.5 GHz fluxes 
vary in the same fashion and that the flux more 
than doubles in the three weeks between the final 
two observations, increasing from 0.23 Jy to 0.47 
Jy at 5 GHz. 



3. The Spectral Index 

The resolution of the VLA at 5 and 8.5 GHz can 
be approximately matched by using the appropri- 
ate weighting schemes, thus allowing for pixel- by- 
pixel measurement of the spectral index a {Sjy cx 
It is important that we understand the ex- 
tent to which the spectral index changes along the 
jet when we calculate the intrinsic brightness (see 
Section |4]) since the Doppler beaming correction 
depends on this quantity. 

To produce spectral index maps we use the uni- 
formly weighted 5 GHz data and the naturally 
weighted 8.5 GHz data. The average beam diam- 
eter (at FWHM) of the uniformly weighted 5 GHz 
images is 370 milliarcseconds (mas), while the av- 
erage beam diameter for the naturally weighted 
8.5 GHz images is 280 mas. The CLEAN mod- 
els for each frequency were convolved with a 350 
mas circular restoring beam so that the two images 
could be compared pixel- by-pixel. The spectral in- 
dex is simply 
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The resulting maps are shown in Figure |H Each 
has been clipped to include data only where the 
estimated spectral index uncertainty is less than 
0.1. In each map, the uncertainty in a ranges be- 
tween 0.05 and 0.1. This estimate includes map 
noise and the 3% calibration uncertainty. The con- 
tours are of total intensity at 5 GHz and increase 
in steps of 2. The kinematic model is displayed 
over each image. The spectral index map for the 
August 24th epoch has not been included. Exces- 
sive errors in the images at this epoch, particularly 
around the core (see the large diagonal streaks in 
Figures [1] and [2]) prevented us from measuring the 
spectral index accurately from these images. 

In addition to matching the image resolutions, 
we must also consider the fact that the UV cov- 
erage between the 5 and 8.5 GHz observations is 
not the same, particularly at short spacings. In 
any inter ferometric observation there is a 'hole' in 
the center of the UV plane that corresponds to the 
minimum distance between antenna pairs. For the 
observations presented here the minimum spacing 
is - 3 kA at 5 GHz and - 6 kA at 8.5 GHz. 

We tested whether missing short spacing flux 
at 8.5 GHz makes a significant change in the spec- 
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tral index images in two ways. First, we sim- 
ply removed visibilities at spacings shorter than 
6 kA from the 5 GHz data to match the minimum 
spacings at 8.5 GHz. Second, we used the AIPS 
task UVSUB to sample a CLEANed C-band image 
with the UV coverage of the X-band observation 
at the same epoch. Each modified UV data set 
was then re-imaged and compared to the original 
C-band image to check for a decrease in intensity 
due to the missing short spacings. In neither case 
were the images or resulting spectral index maps 
significantly changed. 

We find that, within measurement uncertain- 
ties, the spectral index is uniform across the 
source; there is no evidence of significant spatial 
structure or evolution in time. The average spec- 
tral index is found to be 0.68 ± 0.12, 0.73 ± 0.10, 
0.77 ± 0.13, and 0.79 ± 0.15 for June 8th, July 
1st, July 18th, and August 5th respectively. From 
these measurements we conclude that the aver- 
age spectral index is 0.74 ± 0.06, c onsistent with 
the va lue of 0.6 — 0.8 measured by [Stirling et al 



(120041). only slight ly steeper than the 0.6 measured 
by ISeaquist et al.l ^98 0) , and consistent with the 
value of 0.7 used in Paper II. 

We find no evidence for a flattenin g of the spec- 
trum at the core. On VLBI scales, Paragi et al 



(|l999h found an inverted spectrum at the base 
of the jets, which are separated by an AU-scale 
gap that scales inver sely with frequency as pre - 
dicted by the model of lBlandford fc Konigll (|l979h . 
The jets observed in their images are unresolved 
in our images; the entire VLBI scale jet is con- 
tained within the VLA-scale core. Since only the 
base of the miliarcsecond-scale jets has an inverted 
spectrum and the spectrum quickly steepens, we 
should expect the core in our images to have a 
similar spectral index to the rest of the jet. 

There is also no evidence of steepening due to 
radiative losses on this scale. T he half life of a 
synchrotron radiating electron is (jPe Youngll2002L 
eqn 3.57) 



(470 yr) 
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The f iducial field strength is taken from IBcIII 
(|2010l) who measured an equipartition field strength 
in the jets of SS 433 of 10 mG at a distance of ^ \" 
from the core. Therefore the lifetime of the radi- 



ating electrons is much longer than the age of the 
jet material at these scales (less than one year). 
We would have to observe 1100 precessional 
wavelengths along the jet in order to measure any 
steepening at our observation frequencies. This is 
comparable to the time it takes for the jet mate- 
rial to reach the shell of W50, the supernova rem- 
nant i n which SS 433 is contained (lElston fc Baum 
19871 : iDubner et aD [l998). 



4. The Intrinsic Brightness Profiles 

The analysis of radio maps of even mildly rel- 
ativistic jets is complicated by Doppler beaming. 
For SS433, the helical geometry also requires that 
consideration be given to projection effects. For- 
tunately the geometry of SS 433 is well known and 
we are able to correct for these effects in order to 
estimate the intrinsic properties of the jets. In Pa- 
per II we developed methods for calculating the in- 
trinsic brightness and applied them to a very deep 
5 GHz VLA image of SS 433 made from an obser- 
vation in 2003 July. With the data presented here 
we can repeat the same calculations on a series 
of deep images thus allowing us to obtain corre- 
sponding results at several observational epochs. 
Additionally we can study the dependence of the 
intrinsic brightness not only on the age of the ma- 
terial when it emits the radiation that we observe, 
but also on the birth epoch of the material. 

Before describing the corrections we should 
briefly explain what we actually want to measure. 
What does the term "intrinsic brightness" mean 
specifically? If we describe the jets as if they were 
composed of many closely spaced individual com- 
ponents, we would like to know the power emitted 
by each component as measured by a co-moving 
observer. If we describe the jets as continuous, 
we would like to know the emitted power per unit 
length along the jet helix. 

As a result of the (mildly) relativistic jet veloc- 
ity, we must consider the effects of Doppler beam- 
ing. A quantitative measure of this effect requires 
knowledge of the velocity vector at each location 
along the jet helix. Fortunately the necessary ve- 
locities and angles are easily obtained from the 
kinematic model. Additionally, because of the he- 
lical geometry of the jet and the distortion of the 
observed helix due to differences in light travel 
time between different parts of the jet, the amount 



4 



of material that contributes flux to the observing 
beam will differ from location to location. In order 
to determine the contribution per unit length (or 
per component) we must also consider the chang- 
ing projected density on the sky. Again using the 
model, we can calculate the density and normalize 
the measured brightness by this value. 

We will refer to four different time coordi- 
nates that must be defined explicitly. The birth 
epoch, i.e. the date at which a part of the jet, 
as seen by the observer, was emitted from the 
core, will be denoted as 4 (we will always list the 
birth epoch using a Truncated Julian Date, TJD 
= JD - 2454000). The time t denotes the obser- 
vation date for a distant observer (our radio tele- 
scope). We also make use of the "age" of a piece of 
the jet, that is, how long it has been since the piece 
in question emerged from the core. The observed 
age of a piece of the jet is defined as tage = t — tjj 
(assuming a constant velocity). Because the light 
travel time varies along the jet helix, the age of a 
piece of the jet when it emits the observed photons 
is generally not equal to the observed age. The re- 
lationship between the observed age of a piece of 
the jet (tage) and the age of the same jet material 
when it emitted the photons we receive (r) is 

^ _ tgge _ tgge /^x 
~ 1 -/3C0S6> ~ 1-Va:/C ^ ^ 

where is the angle between the local velocity 
vector and the line of sight. This is the same r 
used in Paper II and called "age at time of emission 
of photons". 

If the model information is stored in a list con- 
taining N model jet components equally spaced 
in t^, the total correction factor, Q^n? at the po- 
sition of component i, including both projection 
and relativistic effects, is 

N 

= ^I)-+-e-(^^-)'/(2^'). (4) 

Here D = [7(1 — PcosO)]~^ is the Doppler fac- 
tor, (J is the standard deviation of the Gaussian 
CLEAN beam {FWHM = 2^2 ln2cr), and ARij 
is the angular separation between components i 
and j. We normalize this correction to unity at 
the core. Here n is a parameter that is 2 for a 
continuous jet and 3 if the jet is composed of indi- 
vidual, resolved components. In order to recover 



the "intrinsic brightness profile" of the jet, we sim- 
ply take a reading from the image at the location 
of each model point and divide by this (dimension- 
less) correction factor. 

This correction assumes that each part of the 
jet within a beam area has the same intrinsic lu- 
minosity. This is approximately true for most of 
the jet, but fails at the crossing points of the loops 
on the west. Here, material of very different emis- 
sion ages contributes flux to the same beam. As 
a result, the calculated intrinsic brightness at the 
front of the loop (corresponding to the younger jet 
material within the beam) will be slightly lower 
than it should be, while that of the older compo- 
nents at the rear of the loop will be higher than 
it should be. Without information about how the 
brightness evolves in r and varies with 4 it is im- 
possible to attribute the appropriate amount of 
brightness to each part of the jet. 

The question of the choice of n was addressed 
in some detail in Paper II where we compared pro- 
files from a 2003 VLA image at 5 GHz that had 
been normalized with each value. We determined 
that n = 2 is the appropriate value because the 
jet profiles appear more symmetric and, most con- 
vincingly, the observed west jet profile was very ac- 
curately reconstructed from the intrinsic east jet 
profile. Similarly, here with the 2007 data we find 
that the results are more consistent with a value 
of n = 2 and will use it in the following analysis. 
Figure [5] shows 1/Ci,2, the factor by which each 
sample is multiplied, as a function of r for each 
epoch. 

Figure [6] shows the measured intensity profiles 
from each of the five 5 GHz uniformly weighted 
images. The 8.5 GHz profiles are similar and not 
shown. These have been plotted against emis- 
sion age r rather than the observed age (which in- 
cludes light travel time) because it is r that defines 
the evolution of the system. A Gaussian compo- 
nent having the same dimensions as the restoring 
beam has been subtracted from the images at the 
core prior to measuring the jet brightness which 
is why the profiles fall precipitously to zero as r 
approaches zero. 

The blue and red shaded areas on the plots 
represent the uncertainties of the measurement of 
flux density in the east and west jets, respectively. 
To estimate these uncertainties we consider the 
noise level of the image, a three percent overall 
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flux calibration uncertainty (which is a conserva- 
tive estimate), and the uncertainty in model posi- 
tion. To calculate this last source of uncertainty 
we measure the flux density of the map using 
three models, one created using the mean speed of 
/3 = 0.2647 and the other two with speeds differ- 



i ng by ±10% from the mean. Blundell fc Bowlerl 
([2004) report a 10% fluctuation in velocity, so this 
should give a reasonable measurement at the ex- 
tremes of the possible model positions. We then 
measure the average deviation in the measured 
values from each of these models and take it to 
be the magnitude of uncertainty in the flux density 
due to uncertainty in the position of the jet compo- 
nents. This source of uncertainty dominates over 
the first two; the RMS difference between the pro- 
files measured using different models is ~ 0.5 mJy 
which is roughly an order of magnitude greater 
than the map noise and calibration uncertainty. 

When studying these plots it is helpful to note 
the values of r for which the local velocity vector 
is in the plane of the sky. Table |4] lists these times 
as well as the times when the jet velocity vector 
makes maximum and minimum angles to the line 
of sight. Careful inspection of the observed jet 
profiles (Figure [6]) shows that at the sky crossing 
times the brightness of the east and west jets is 
very similar, as expected since the Doppler factors 
and projection at those times are the same. At 
other times there appear to be asymmetries in the 
measured brightness, and the east and west pro- 
files in general have rather different shapes. For 
instance, there is a large difference between the 
east and west intensity values between 50 and 100 
days in the early observations. The asymmetries 
and features in the curve, such as the sharp break 
at 100 days in the June 8th profile, appear in the 
profiles of each epoch at later times as the sum- 
mer goes on. The aforementioned difference, for 
instance, moves outward to 100 to 170 days by 
August 24th. 

Figure [71 shows the intrinsic brightness profiles. 
The corrected profiles appear much more symmet- 
ric, with the only major difference being a small 
dip followed by a large bump in the west jet pro- 
file. This is seen at r ~ 150 — 170 days in the 
June 8th profile and appears at later times in sub- 
sequent images. As described above, this apparent 
mismatch results from components of very differ- 
ent ages contributing fiux to the map at the loop 



crossing point in the west jet (see, for instance, the 
region at relative right ascension —1.2'^, declina- 
tion 0.5'^ on the July 1st image in Figure [1]). The 
minimum of the small dip coincides with the front 
of the jet cone while the maximum corresponds to 
the value measured at the location of the "older" 
model point at the rear of the jet cone. 

There are other, smaller differences between the 
east and west jet intensities in the corrected pro- 
files. For instance, around r = 100 days in the 
July 1st profile, the intrinsic brightness measured 
in the west jet is 20% higher than that measured 
in the east jet. This difference is slightly outside 
of the estimated uncertainties. However, for this 
value of r there is a difference in tb between the 
east and west jet of ~ 10 days. There are similar 
differences at the other epochs as well, but they 
all appear when there is a significant difference in 
birth epoch between the east and west jet compo- 
nents at a given r. In each profile, where At 5 = 
the profiles match quite well. We conclude that 
the apparent difference between profiles is indica- 
tive of intrinsic fiux variation with ^5. We investi- 
gate this further in the next sections. 

4.1. Evolution of the Intrinsic Brightness 

When examining the intrinsic jet profile of a 
single epoch as we did in Paper II it is immedi- 
ately clear that the curve is not well described by 
a simple mathematical function. From the 2003 
July 11 VLA observation we found three regions, 
each best described by different power law or expo- 
nential functions. In the range 50 < r < 150 days, 
the intrinsic brightness curve was adequately fit 
by a power law with exponent —1.8 or an expo- 
nential function with a half-life of ~ 40 days. The 
curve then fiattened, displaying a roughly constant 
intrinsic brightness from 150 < r < 300 days. Be- 
yond ^ 300 days the intrinsic brightness fell off 
more quickly, and was best fit by a power law with 
an exponent of roughly —4 or an exponential func- 
tion with a half-life of ~ 80 days. 

Hjellming fc Johnston I (j 19881 ) also noted a 



break in the jet brightness profile, with the rate 
of decay increasing sharply at tage = 100 days. 
Based on their model of a precessing conical jet 
(discussed below), they attributed this behavior 
to a change in the transverse jet expansion rate. 
They argued that the jets initially undergo con- 
fined expansion and then, after about 100 days. 
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the jets expand freely. 

Similar breaks are apparent in the profiles of the 
2007 data presented here. From our multi-epoch 
data set we see that these features propagate out- 
ward from epoch to epoch. The brightness curves 
from the 2007 images are quite different from that 
of the 2003 epoch. Figure [8] compares the intrin- 
sic east jet profiles of the 2003 July 11 epoch (Pa- 
per II) and the 2007 July 18 epoch. On these dates 
the source was at a similar precessional phase, 0.92 
in 2003 and 0.96 in 2007. 

Observations of the jet profile from a single 
epoch are insufficient for characterizing the evo- 
lution of the jet brightness since, in general, the 
brightness of any piece of the jet might be ex- 
pected to depend on two time parameters. These 
are ^5, the birth epoch, and r, the age of the com- 
ponent when it emits the photons that we observe. 
The dependence on tb reflects the varying level of 
power injected into the jets at the core, and the 
dependence on r reflects the evolution of individ- 
ual pieces of jet. From the jet profile at a single 
epoch we can measure the brightness of the jet for 
many different values of 4, but we only measure 
the brightness for one value of r at each value of 
ti). With the multi-epoch observations presented 
here, however, we have a measure of the intrinsic 
brightness at up to five values of r for each 4 (ten 
if we assume east /west symmetry). 

We will assume that the functional form for the 
brightness at each location along the jet can be 
factored into parts that each depend only on one 
of these times, 

I{h,T)=C{h)f{T). (5) 

Under this assumption we can measure /(r) if we 
simply measure the intrinsic brightness at a point 
with a given at each epoch. In Figure [9] we show 
the evolution of the brightness, measured from the 
east jet, for two different values of ^5. The bright- 
ness values in Figure [9^ are measured at tb = 143 
days, which corresponds to r = 111 days and the 
image location (0.9, —0.2)'' on June 8th. In Fig- 
ure [9b the brightnesses are measured at = 179 
days, which corresponds to r = 83 days and the 
image location (0.5,-0.1)'' on June 8th. We have 
chosen points that are well separated on the image 
and whose intensities have a relatively low uncer- 
tainty. Power law (r~*) and exponential (e"'^/'^ ) 
fits are shown as dashed red lines and solid black 



lines, respectively. For the values measured when 
h = 143 days we find 5 = 2.4 ± 0.5 and r' = 58 ± 6 
days. For U = 179 we find s = 1.9 ± 0.8 and 
r' = 58 ±10 days. Both models appear to describe 
the data equally well. The results are consistent 
with a constant fitting parameter for either model, 
although the prediction bounds are quite large due 
to the limited number of data points used in the 
fit. 

We can improve our description of the evolu- 
tionary behavior if we repeat the measurement 
above for all values of 4. We assume, as suggested 
by Figure [HI that /(r) can be represented by ei- 
ther a power law or exponential function. One of 
the intensity measurements at a given will be 
denoted as /o(^6) and the corresponding emission 
age as tq . Then it follows from equation [5] that 
the intensity and emission age at other times are 
related to these by either 

if /(r) is a power law, or 

— = e-^^-^°)/^' (7) 
^0 

if /(r) is an exponential. These ratios have the 
benefit of being independent of tb and therefore 
free from effects of variation in the core. With 
these ratios we can include all measured inten- 
sity values in order to better constrain our best 
fit parameters, s and r', and possibly distinguish 
between the two models. 

Figure [To] shows the ratio /(t^, r)//o(t5, tq) as a 
function of t/tq. This is plotted on a log- log scale 
so a power law will appear linear. This includes in- 
trinsic brightness values from both the east (shown 
in black) and west (shown in red) profiles. Mea- 
surements from the west jet in the loop crossing 
region, and from regions where the core accounts 
for more than 1% of the total flux, have been dis- 
carded. The best fit of equation [6] to the east jet 
data, shown as a solid line in the figure, yields an 
exponent of s = 2.0 ±0.1. The best fit to the 
west jet data, shown as a dashed line in the figure, 
yields an exponent of s = 1.8 ± 0.1. The (co- 
efficient of determination) values for the fits are 
0.76 and 0.70 for east and west, respectively. 

Figure [TT] shows the ratio 7(^5, r)//o(t5, tq) as a 
function of r — tq. On this plot an exponential will 
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appear linear. The data have been selected using 
the same criteria as in Figure [TOl The best fit of 
equation [71 to the east jet, shown as a solid line in 
the figure, gives t' = 57.6 ±1.4 days. The best 
fit to the west jet, shown as a dashed line in the 
figure, gives t' = 54.3 ± 1.9 days. The exponential 
fits are better than the power law fits for both 
jets; the values for the fits are 0.93 and 0.90 for 
east and west, respectively. The fits are consistent 
with the east and west jets having the same value 
of r', indicating that the external environments 
surrounding each jet are comparable. 

We note that there is some vertical scatter in 
each of these figures, and especially in Figure [TOl 
This is not due to noise, rather it is the result of 
one of several possible systematic effects. For ex- 
ample, if the assumed functional forms for /(r) are 
not valid then this kind of scatter will be expected. 
In Figure [To] we suspect that a large portion of the 
scatter is due to the inadequacy of the power law 
model. Another possibility is that this is due to 
a breakdown in the assumption made in equation 
[SI e.g. if To varied slightly with t^. Nevertheless, 
from these figures we find that these kinds of sys- 
tematic effects are secondary to the general, simple 
trends analyzed above. 

We conclude that the evolution of jet bright- 
ness in r, within the range 50 < r < 300 days for 
both the east and west jets, is well described by 
an exponential function with t' = 55.9 ±1.7 days. 
This is a significantly longer exponential time con- 
stant than the value of 13 — 21 days reported by 
jjowett fc Spencer! ( 1995 ). We note, however, that 
they were measuring the brightness evolution as a 
function of tage rather than r and thus have not 
taken into account the varying amounts of light 
travel time from component to component. 

4.2. History of the Core Flux Variability 

If /(r) has been determined then we can recover 
the proportionality constant C for a given birth 
epoch simply: 



C{h) = I{h,T)/f{T). 



(8) 



Here we assume /(r) = e~^l^ and we use t' — 
55.9 days. Figure [12] shows this brightness coef- 
ficient as a function of birth epoch for both the 
east (shown in black) and west (shown in red) jets. 
This includes data from all five epochs. The val- 
ues measured for the west jet cover a much smaller 



range in because we have to exclude data mea- 
sured from the "loop". 

The coefficient increases by more than a factor 
of five over the course of 50 days between TJD 100 
and 150. Where we have measured C in both jets, 
the results are consistent with the east and west 
jet values being the same. 

The general trend is consistent between 5 and 
8.5 GHz as well. In Figure [13] we show the ratio of 
the east jet coefficients measured from the 5 GHz 
profiles to those measured from the 8.5 GHz pro- 
files {C^GHzith) / C^,^GHz(th))- The ratio remains 
roughly constant throughout the range in ^5. Also 
note that, after accounting for the different beam 
dimensions, this ratio is consistent with the core 
flux ratio shown in Figure [H The agreement be- 
tween 5 and 8.5 GHz coefficients indicates vari- 
ability with constant spectral index. 

The measurement of C{ti)) provides a history 
of the relative power injected in the jet and could 
be used as a post-diction of the history of the ra- 
dio core variability over the year preceding the ob- 
servations. This interpretation suggests that the 
core fiux would have increased dramatically over 
the course of a few weeks. This kind of outburst 
is not unprecedented, with a similar fiare being 
detected during a series of VLA observations of 
SS433 in the summer of 2003 (Roberts et al., in 
preparation). We are not aware of any monitor- 
ing campaigns that span the period from mid 2006 
to mid 2007, so this post-diction is not verifiable. 
Application of this procedure to VLA data taken 
in collaboration with a radio monitoring campaign 
would be required for confirmation. 

5. Discussion 

5.1. Comparing the Brightness Evolution 
to Radio Source Models 

Now that we have a measurement of the bright- 
ness decay rate, we can compare the observed rate 
with model predictions. The simplest model is 
that of a spherically expanding radio source pro- 
posed by Ivan der LaanI (|l966^ . In the case of 
SS 433 we can think of the je ts as a series of spher- 
ical blobs dVermeulen et al ][i987). In the van del 



Laan model, Sj^ (x r~^^, where p = 1 + 2a is the 
exponent of the assumed electron energy distri- 
bution and r is the radius of the source. In the 
case of ballistic expansion, where r oc r, and since 
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here p 2.4, this model predicts that the bright- 
ness will decay as a power law with an exponent 
of —4.8, which is much faster than we measure. 
If instead the jet expands at a slower rate, with 
r cx r^/^, then the decay would be proportional to 
r~^-^ which is comparable to the power law fits to 
the data. An exponential decay of the brightness 
in time would require that the radius of the blobs 
expands exponentially. If the brightness decays 
with an exponential time constant of 55.9 days 
and p = 2A then the radius will increase with an 
exponential time constant of 268 days. 

A mo del of a precessin^ conical outfl ow was pre- 
sented in'Hjellm ing fc Johnstonl ( 1988 ^ to describe 
the observed jet proflle in SS 433 and other pre- 
cessing jet systems. In the optically thin limit they 
flnd that the contribution to the flux density from 
each segment of the jet, dz, is dSu{z) cx z~^dz, 
where z is the distance away from the core (along 
the local velocity vector) and m = {7p—l)/{6-\-6S). 
The parameter (5 is 1 if the jet expansion is con- 
flned or if it expands freely. Since the speed 
of the jet is assumed to be constant, z cx r, so 
this also provides a relationship between the flux 
per unit length and time. For p 2.4, this re- 
sults in Si, cx r~^"^ or r~^-^ for the free and con- 
flned expansions, respectively. Comparison with 
our power law flts would suggest that the jets ex- 
pand at a rate between these, with the east jet 
expanding more freely than the west. However, 
we do not believe that this model appropriately 
describes the jets of SS433 because of an appar- 
ent flaw in the assumed geometry of the model. 

One of the key assumptions of themodel pro- 
posed by iHjellming fc Johnstonl ( 1988[ ) is that 
there are no gradients in the velocity along the 
jet locus. Therefore the jets do not expand longi- 
tudinally but only laterally (2-D expansion). This 
assumption makes sense for a straight jet. For a 
helical jet, however, even though the magnitude of 
velocity is nearly the same for all parts of the jet 
and does not change with distance from the core 
as best we can measure (see Section [5^ . the di- 
rection of the velocity vector varies from location 
to location. As a result, one component of the jet 
does not follow its nearest neighbor. Instead their 
trajectories diverge and expansion along the jet 
helix is inevitable. 



5.2. Comparing the Images to the Kine- 
matic Model 

We know that the jet material is traveling bal- 
listically through an external medium that is com- 
posed of the disk wind and other material within 
the surrounding supernova remnant. The jets will 
interact with the external environment and decel- 
erate due to the resulting ram pressure. The jets 
observed in our radio images however, do not ap- 
pear to have decelerated measurably. This is made 
clear by comparing the images to the kinematic 
model which successfully predicts the location of 
the jet out to 4'' from the core. The parameters 
that we use to generate this model (apart from 
the distance to the source) have been derived in- 
dependently from optical emission line data. The 
emission lin es originate from a region sm aller than 
- 10^^ cm ([Davidson fc McCravl [T98Qh while the 
arcsecond scale radio maps reveal the jets at a 
scale of 10^^ cm. The model curve is computed 
assuming constant velocity. 

(|2004l ) 



Stirling et al 



found evidence for decel- 
eration of 0.02c per period by comparing the kine- 
matic model locus with radio images made using 
VLA and MERLIN data. It should be noted, how- 
ever, that they assume a distance of 4.8 kpc where 
we adopt the more recently verifled distance of 5.5 
kpc from iBlundeU fc Bowlerl |2004l). Calculating 
the model positions using a shorter distance will of 
course result in a larger projected size of the jets 
and accounts for the discrepancy between model 
and image that they noted. Inclusion of a 0.02c 
per period deceleration in the model calculation 
results in a compressed appearance of the jets that 
we do not observe. Moreover, the kinematic model 
computed with a constant velocity flts our images 
equally well throughout the summer of 2007. If 
the jets were decelerating we would expect the 
model to systematically lead the observed helix 
at increasing angular distance from the core. 

The excellent agreement between the radio im- 
ages and model curves is indicative of the stabil- 
ity of the velocity across a broad range of length 
scales. We can use the agreement between the 
kinematic model and our radio images to put lim- 
its on the deceleration of jet material and use this 
to estimate the ratio of mass density between the 
jet and surrounding medium. For simplicity, we 
assume that the jet is composed of spherical blobs 
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of plasma. Each blob has radius r. If we relate 
the ram pressure acting on the blobs to the accel- 
eration we find the following relationship between 
the density of the internal and external media, 



Pi 

Pe 



(9) 



where pe is the external density, pj is the jet den- 
sity, V is the blob velocity relative to the external 
mediumjl and is a shape factor of order unity 
(we set = 1). The density pj will certainly de- 
crease as the jets cool and expand. The density 
of the surrounding medium, will also likely de- 
crease with increased distance from the binary sys- 
tem. For simplicity we assume that the densities 
decrease in such a way that their ratio is constant. 

Any acceleration would result in a displacement 
from the un-accelerated position (represented by 
the model curve) at a given time. The maximum 
displacement is Ad < \vt^ • We will assume that 
this displacement projects to an angular distance 
that is less than the spatial resolution of our im- 
ages since otherwise we would be able to detect it. 
An image resolution of 0.3 arcseconds corresponds 
to a linear size of 2 x 10^^ cm, assuming a distance 
to SS433 of 5.5 kpc. Substituting the maximum 
displacement for the acceleration we find 



Pe 8rAd ' 



(10) 



We measure the jets clearly out to an age of ~ 350 
days. For a ballistically expanding jet with a cone 
half-opening angle of 0.6° (Marshall et al. 2002j) 
this ratio is pj / pe ^ 300. For a slower expansion 
rate the ratio will be larger. 

6. Conclusions 

We have presented analysis of a sequence of five 
deep observations of SS 433 made over the summer 
of 2007 using the VLA at 5 and 8 GHz. The main 
results are as follows. 



^Begelman et al] ([lOSOT ) find that the disk wind travels with 
a velocity of 1500 km/s, which is only ~ 2% of the jet 
speed (~ 8 X 10^ km/s). Thus we assume that the external 
environment is stationary relative to the jets. 



We measure the spectral index to be 0.74 ± 
0.06, consistent with previously reported 
values, and find no significant variations ei- 
ther across the source or with precession 
phase. 

All the results of Paper II (based on a single 
very deep 5 GHz observation made in 2003) 
are confirmed, but with the sequence of five 
closely spaced observations, the possibilities 
for analysis are greatly increased. 

The jet intrinsic brightness profiles in 2007 
are significantly different from those in 2003, 
even though the precession phases are sim- 
ilar. We attribute this to variability at the 
core, injecting varying power at the base of 
the jets. 

As in Paper II, in every image the profiles of 
the east and west jets are remarkably similar 
if projection and Doppler beaming are taken 
into account. 

The sequence of five images allows us to dis- 
entangle the evolution of individual pieces of 
jet from variations of jet power originating at 
the core. We find that the brightness of in- 
dividual pieces of jet fades as an exponential 
function of age (or distance from the core) 
rather than a power law. The exponential 
time constant is 55.9 ± 1.7 days. Our results 
are consistent with the brightness of the two 
jets evolving in the same way. 

There is also significant variation (by a fac- 
tor of at least five) in jet power with birth 
epoch, with the east and west jets varying in 
synchrony. The variation is also consistent 
between 5 and 8.5 GHz. 

The lack of deceleration between the scale of 
the optical Balmer line emission (10^^ cm) 
and that of the radio emission (10^^ cm) re- 
quires that the jet material is much denser 
than its surroundings. We find that the den- 
sity ratio must exceed 300:1. 
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Table 1 
Observational Summary 







# Antennas (VLA/EVLA) 




Date 


TJD^ 


C-Band 


X-Band 


Precessional Phase^ 


2008 June 6 


260 


16/4 


16/9 


0.72 


2007 July 1 


283 


17/5 


17/8 


0.86 


2007 July 18 


300 


16/6 


15/9 


0.96 


2007 Aug. 5 


318 


14/7 


16/9 


0.07 


2007 Aug. 24 


337 


15/5 


16/7 


0.19 



^Phase is zer o whe n the east jet is maximally redshifted following 
Vermeulen et al.l (|l993h . 



^Truncated Julian Date in days, TJD = JD - 2454000 



Table 2 
C-BAND Image Details 



Date 


BMAJ 


BMIN 


BPA 


^max 




RMS 




arcsec 


arcsec 


Degrees 


mJy/Beam 


mJy/Beam 


mJy/Beam 


2007 June 8 


0.404 


0.340 


-68 


329 


0.12 


0.04 


2007 July 1 


0.381 


0.335 


-22 


359 


0.10 


0.03 


2007 July 18 


0.410 


0.348 


-69 


261 


0.11 


0.04 


2007 Aug. 5 


0.414 


0.345 


-88 


225 


0.14 


0.05 


2007 Aug. 24 


0.356 


0.330 


35 


467 


0.19 


0.05 



Table 3 
X-BAND Image Details 



Date 


BMAJ 


BMIN 


BPA 


^max 




RMS 




arcsec 


arcsec 


Degrees 


mJy/Beam 


mJy/Beam 


mJy/Beam 


2007 June 8 


0.219 


0.208 


53 


211 


0.08 


0.03 


2007 July 1 


0.219 


0.201 


-31 


221 


0.08 


0.03 


2007 July 18 


0.224 


0.209 


63 


141 


0.09 


0.03 


2007 Aug. 5 


0.238 


0.205 


-73 


129 


0.07 


0.02 


2007 Aug. 24 


0.240 


0.213 


-49 


361 


0.11 


0.04 
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Table 4 

Profile Plot Reference Times (t in days) 



Obs. Date 


ee = o^ = 90° 


6e = min. 


Oe = max. 


6^ = min. 


Oy^ = max. 


2007 June 8 


91, 142, 253, 304 


41, 231, 420 


110, 267 


121, 291 


30, 173, 315 


2007 July 1 


114, 165, 276, 327 


67, 257, 447 


133, 289 


146, 315 


52, 192, 335 


2007 July 18 


20, 131, 182, 293, 344 


88, 278 


150, 305 


163, 333 


66, 208, 350 


2007 Aug. 5 


38, 149, 200, 311, 362 


107,298 


12, 168, 333 


12, 182, 352 


81, 223 


2207 Aug. 24 


57, 168, 219, 330, 380 


131, 320 


30, 185, 341 


32, 202, 370 


98, 240 
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Fig. 1. — C-band uniformly weighted images of total intensity. Contours are in steps of a/2. The restoring 
beam is represented by a cross in the lower right of each image. Overlayed on each figure is the kinematic 
model with blue and red line segments representing the retreating and oncoming jet material, respectively. 
Detailed information about contour levels, map noise, etc. is found in Table El (a) 2007 June 8, (b) 2007 
July 1, (c) 2007 July 18, (d) 2007 August 5, (e) 2007 August 24. 
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Relative Right Ascension (arcsec) 




Fig. 2. — X-band uniformly weighted images of total intensity. Contours are in steps of a/2. The restoring 
beam is represented by a cross in the lower right of each image. Overlayed on each figure is the kinematic 
model with blue and red line segments representing the retreating and oncoming jet material, respectively. 
Detailed information about contour levels, map noise, etc. is found in Table [3l (a) 2007 June 8, (b) 2007 
July , (c) 2007 July 18, (d) 2007 August 5, (e) 2007 August 24. 
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Fig. 3. — The core flux as a function of TJD for (black) 5 GHz and (red) 8.5 GHz. The error bars represent 
a 3% calibration uncertainty. The 5 and 8.5 GHz fluxes vary in the same manner, and each roughly doubles 
between 2007 August 5 and 2007 August 24. 
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Fig. 4. — Maps of the spectral index a (S^y ex The contours are of C-band total intensity with steps 

of 2 between levels. The restoring beam was specified to be circular with full- width half-max 0.350'^ The 
kinematic model is displayed over each image, (a) 2007 June 8, (b) 2007 July 1, (c) 2007 July 18, (d) 2007 
August 5. 
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Fig. 5. — The multiplicative correction factors. The intrinsic brightness is the product of this correction 
factor and the observed total intensities, (a) 2007 June 8, (b) 2007 July 1, (c) 2007 July 18, (d) 2007 August 
5, (e) 2007 August 24. 
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Fig. 6. — Observed total intensity profiles from the C-Band observations as a function of the emission age. 
The colored regions indicate the estimated uncertainties; blue for the east jet and red for the west. The 
core was subtracted from each image prior to measuring the profiles. The bump in the west jet profiles that 
moves from '^150 to 250 days is where the projected helix crosses itself (see text), (a) 2007 June 8, (b) 2007 
July 1, (c) 2007 July 18, (d) 2007 August 5, (e) 2007 August 24. 



20 






Fig. 7. — The intrinsic brightness profles of the C-Band observations as a function of the emission age. The 
colored regions indicate the estimated uncertainties; blue for the east jet and red for the west. The core was 
subtracted from each image prior to measuring the profiles. The bump in the west jet profiles that moves 
from ^150 to 250 days is where the projected helix crosses itself (see text), (a) 2007 June 8, (b) 2007 July 
1, (c) 2007 July 18, (d) 2007 August 5, (e) 2007 August 24. 
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Fig. 8. — Comparison between the east jet intrinsic brightness profiles of 2003 July 11 and 2007 July 18. 




Fig. 9. — The evolution of the C-band intrinsic brightness for two pieces of the eastern jet emitted at different 
birth epochs (^5). The uncertainty is estimated using the procedure employed for Figures [6] and The best 
power law and exponential fits are shown as dashed red and solid black lines, respectively, (a) Fits for 
h = 143: = 58 ± 6 days and s = 2.4 zb 0.5, (b) fits for tt = 179: = 58 ± 10 days and s = 1.9 ± 0.8. 
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Fig. 10. — The normalized intrinsic brightness, I{ti},r)/Io{tb^ro), as a function of t/tq. The trend wih be 
Unear if the brightness is a power law in r. Values measured from the east and west jet are shown as black 
dots and red x-marks, respectively. The best power law fit is shown as the solid line for the east jet and the 
dashed line for the west jet. For the east jet s = 2.0 ±0.1, and for the west jet s = 1.6 ± 0.1. 
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Fig. 11. — The natural logarithm of the normalized intrinsic brightness, 7(^5, r)//o(t5, tq), as a function of 
the difference (r — tq). The trend will be linear if the brightness depends exponentially on r. The values 
measured from the east and west jet are shown as black dots and red x-marks, respectively. The best 
exponential fit is shown as the solid line for the east jet and the dashed line for the west jet. For the east 
jet t' = 57.6 ±1.4 days, and for the west jet t' = 54.3 ±1.9 days. 
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Fig. 12. — The birth epoch dependent part of the intrinsic brightness, C{tb)^ as a function of birth epoch mea- 



sured in TJD. We compute C assuming that the r dependent part is an exponential, I{tb^r) = C{th)e 



-T I t' 



with t' = 55.9 days. The values derived for the east and west jets are shown in black and red, respectively. 
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Fig. 13. — The ratio of C{tb) as measured from the east jet in the 5 and 8.5 GHz images. The shaded 
indicates the uncertainty. 
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